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ABSTRACT 

Wc analyse for the first time the radial abundance gradients of the disc stars of a 
disc galaxy simulated with our three dimensional, fully cosmological chemodynamical 
galaxy evolution code GCD+. We study how [Fe/H], [N/0], [0/Fe], [Mg/Fe] and [Si/Fe] 
vary with galactocentric radius. For the young stars of the disc, we found a negative 
slope for [Fe/H] and [N/0] but a positive [0/Fe], [Mg/Fe] and [Si/Fe] slope with 
radius. By analysing the star formation rate (SFR) at different radii, we found that 
the simulated disc contains a greater fraction of young stars in the outer regions, while 
the old stars tend to be concentrated in the inner parts of the disc. This can explain 
the positive [a/Fe] gradient as well as the negative [N/0] gradient with radius. This 
radial trend is a natural outcome of an inside-out formation of the disc, regardless of 
its size and can thus explain the recently observed positive [a/Fe] gradients in the 
Milky Way disc open clusters. 

Key words: Galaxy: disc — Galaxy: kinematics and dynamics — galaxies: Interac- 
tions — galaxies: Formation — galaxies: evolution — galaxies: abundances 



1 INTRODUCTION 

In recent years, there have been extensive efforts devoted to 
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these studies has been to find out how our Galactic 
disc formed and evolved with time. Around the solar 
neighborhood, good progress has been made recently, 
including the derivation of age information for sev- 



eral clusters fe.e. Pancino et al. 


2010|) and individual 


stars fe.e. Nordstrom et al. 
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iJoreensen & LindeErenI l2005l: iHolmberE et al.ll2007l. l2009f). 



However, further out, at larger galactocentric radii, there 
is a striking lack of any large high quality datasets. To 
unravel the formation history of the Galactic disc, we need 
to know the abundance trends in the stars along the disc 
complemented with accurate age information. 

In the near future ESA's Gaia mission should provide 
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exceptional positional and proper motion information for 
up to one billion stars in the Milky Way. The vast majority 
(80%) of these stars will lie in the disc. For bright stars, the 
Gaia RVS will determine the abundances of iron and the 
Q-elements. In addition, future and proposed ground-based 
projects including the Apache Point Observatory Galac- 
tic E volution Experiment, APOGEE (jAUende Prieto et al.l 
l2008f ) and the High Resolutio n Multi-object Echelle Spec- 
trogr aph (HERMES) survey (|Freeman fc Bland-HawthornI 
fiooi) will provide more detailed and accurate chemical 
abundance determinations which will complement the Gaia 
data. 

The chemical properties of stars are important to under- 
stand the formation history of the disc. Chemical elements 
heavier than boron are the end products of stellar evolution. 
The so-called a-elements and iron (Fe) are of particular in- 
terest since it is known that they are produced primarily 
in Type II (SNe II) and Type la (SNe la) supernovae re- 
spectively. SNe la and SNe II have different timescales and 
thus studying the abundance ratios of the a-elements with 
respect to Fe gives unique fossil information on the past con- 
ditions and evolution of the galaxy. In addition to studying 
various [a/Fe] ratios, we also look at how the ratio of [N/O] 
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varies with galactocentric radius as this can also be used as 
a cosmic clock. Good cosmic clocks are obtained by taking 
the ratios of elements produced by different stellar masses 
and therefore on different time-scales. Nitrogen is primarily 
produced by intermediate mass stars with longer lifetimes 
whereas oxygen is produced by SNe II progenitor stars which 
have a shorter lifetime. In any case, the actual abundance 
gradients and ratios of different elements will vary depend- 
ing on the chemical evolution histories at different radii and 
therefore on the building-up history of the disc. 

The chemical evolution at different radii of the disc is 
studied theoretically using pu r e chemical evolution mod - 



els (e.g. IChiappini et all Il997l . I2001I : iRenda et al.l l2005bh . 



IChiappini et al.l ( 19971 ) suggested that the Milky Way forms 
primarily out of two infall episodes; the first giving rise to 
the halo and subsequently the bulge, and the second pro- 
ducing the disc via a much slower infall of primordial gas 
preferentially accumulating faster in the inner compared to 
the outer regions of the disc. This scenario is known as the 
"inside-out" mechanism for disc formation. It is one of the 
possible mechanisms which can reproduce and explain the 
observed properties of the Milky Way. However in order to 
be able to draw firmer conclusions, it is necessary to have 
obse rvational data for the outer regions of the Galactic disc 
ifChla ppini ct al. 200l]). 

Inside-out formation may lead to some specific evolution 
of the abundance gradients. How do the abundance gradi- 
ents evolve within the disc? There is still a lack of agreement 
between di fferent authors; for so me the gradients steepen 
with time ijChiappini et al.lll997^■ whilst t o others the gra- 
dients flatten with time |Hou et al.l I2OO0I : iDaflon fc Cunhal 
I2OO4I : iMaciel fc Costal l2009h. Note, however, that some au- 
thors have recently suggested that the traditional chemica l 
evolution models such as those of IChiappinl et all (I2OOII ') 
may need to be revised as they do not sufliciently consider 
the irn portant effects of radial mixing ijSchonrich &: BinnevI 
l2009al ). 

In recent years, there has been progress in modelling 
disc galaxy formation in a Cold Dark Matter (CDM) 
Universe, u sing three dimensional numerical simulations 



llKatzlll992l:jSteinmetz fc Mulleijll994i:,Bekki fc Chiba,200a: 
Abadi et al.ll2003allbl: iBrook et al.1 l2004bl: iGovernato et al l 



20071 : lOkamoto et al.l [20081 : IScannapieco et al.ll2008l . 120091 ). 

Some studies include both SNe II and SNe la and discuss the 
details of the chemical pro perties in the simul a ted galaxies in 
an isolated halo collapse (iRaiteri et al.l 19961: lBercziklll999l) 

' 200a: ' 



or hie r archical clustering (iBrook et al 



2005al : IScannapieco et ahl 
20081 ) scenario. However, to date, no three dimensional fully 
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20051 : iMartmez- Serrano et al 



cosmological chemodynamical simulations have looked at ra- 
dial abundance trends with age of disc stars. This is the first 
study to do so and it will be very interesting to see how our 
results will compare with the literature, both with theory 
and observations. 

In this paper, we analyse the chemodynamical prop- 
erties of the disc stars in a Milky Way size disc galaxy in 
our ACDM cosmological simulation. Unfortunately, the sim- 
ulated galaxy is not a late-type disc galaxy, such as the Milky 
Way. Nevertheless, since the detailed chemical distribution 
of disc stars is only available to us in the Milky Way, we com- 
pare our simulation with the Galaxy, and discuss what we 



could learn about the formation and evolution of a general 
disc component. 

The outline for this paper is as follows. In Section 2 we 
describe our numerical simulation and define our disc stars. 
In Section 3 we present the results of our chemodynamical 
analyses. We initially study the properties and differences 
between accreted and in situ stars within the disc in Section 
3.1. In Section 3.2 we analyse any radial abundance trends 
in the disc and its relation to the age of the disc stars with 
the aim of unravelling the formation mechanisms of our disc. 
Finally we present our conclusions in Section 4. 



2 THE CODE AND MODEL 

To simulate our galaxy, we use the original galac- 
tic chemodynamic a l evo lution code GCD+ developed by 
iKawata fc GibsonI (|2003l ). GCD+ is a three-dimensional 
tree iV -body /smooth ed particl e hyd r odynamics code (ILuc^- 



I977I: iGingold fc M onagha n [l977l: iBarnes fc HutI Il986. 
Hernguist fc Katj Il989l : IKatz et al.1 |l996^ that incorpo- 



rates self-gravity, hydrodynamics, radiative cooling, star for- 
mation, supernova feedback, and metal enrichment. GCD+ 
tak es into acco unt chemical enrichm ent by both SNe II 
tWooslev fc Wea ver 1995) and SNe la (jlwamoto et al.lll999l : 
iKobavashi et al.1 l2000i) and mass loss from int ermediate- 



mass stars Ijvan den Hoek fc GroenewegenI [19971 ). and fol- 
lows the chemical enrichment history of both the stel- 
lar and gas components of the system. As described in 
[Kawata fc Gibson ( 2003) GCD+ takes in to account the metal- 
licity dependence of the age of stars l|Kodama fc Arimotd 
Il997f ). the metal dependent yields from SNe II and mass 
loss from intermediate mass stars. 

Radiative cooling, which depend s on the metallicity of 
the gas (derived with MAPPINGSIII: [Sutherland fc Dopital 
Il993l ) is taken into account. The cooling rate for a gas with 
solar metallicity is larger than that for gas of primordial 
composition by more than an order of magnitude. Thus, 
cooling by metals should not be ignored in numerical sim- 
ulations of gala xy formation (.Kallander fc Hultman 19981 : 
iKav et al]|2000l) . However, we ignore the heating effect of 
the cosmic UV background radiation and UV radiation from 
hot stars for simplicity. 

Star formatio n is m odelled usi ng a metho d similar to 
that suggested bv iKatj (,1992) and iKatz et al.1 (il996 ). For 
star formation to occur, the following three criteria must 
be satisfied: (i) the gas density is greater than some crit- 
ical density; (ii) the gas velocity field is convergent and 
(iii) the Jeans instability condition is satisfied. Our SFR 
formula corresponds to the Schmidt law. We assume that 
stars are distributed according to the lSalpeterl ll 19551) initial 
mass function (IMF) . For more details see Kawata fc GibsonI 
(2003). Note that we assume only thermal energy feed- 
back from SNe. It is known that the thermal energy ieeA- 
back mode l has a negligible eff ect on galaxy formation (e.g. 
lKatj|l992l : iBrook et aLll2004al ). In real galaxies, the effects 
of SNe feedback are likely stronger, and could affect the 
chemical evolution in th e disc component, as discussed in 
IScannapieco et ahl (|2008l ). However, there is still no model 
that can convincingly explain how SNe feedback affects 
galaxy formation and evolution. In this paper rather than 
exploring a variety of feedback modelling, we have imple- 
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Table 1. Simulation parameters 
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Figure 1 . Histogram showing the rotational velocities of stars in 
our disc region. We defined disc stars to have rotational velocity 
between 150-350 kms~^ . The circular rotation speed in the central 
regions of our disc, at a radius of 7 kpc was 260 kms~^. 



mented one simple feedback model and study how chemical 
evolution takes place under this assumption. Our simulation, 
therefore, do es not include any outflow, such as a Galac- 
tic fountain 11 Bregma ril 119801 ). Although some authors (e.g. 
iRecchi et al.|[200ll ) suggest that SNe la dump more energy 
to the surrounding ISM than SNe II, we assume the same 
energy per supernova for SNe II and SNe la for simplicity. 

The galaxy si mulated here is from the sample of 
iRahimi et al] (I2OI0I '). referred to as "Gall". Gall is a high 



resolution version o f galaxy "Dl" in iKawata et al.l l|2004 
iRahimi et al.l (|201Clil analyse a second galaxy ("Gal2"), but 
Gal2 has a lower resolution and a less prominent disc, and 
thus we focus only on Gall in this paper. We used the multi- 
resolution technique in order to maximise the mass resolu- 
tion wi thin the regions whe re the disc progenitors form and 
evolve (|Kawata et al.l [20041 '). Here, only the high resolution 
region includes the gas particles, and therefore star forma- 
tion. 

We summarise th e properties of our si mulated galaxy in 
Table 1 adapted from IRahimi et al] (|2010l ). The second col- 
umn represents the virial mass; the third column, the virial 
radius; Columns 4 and 5 represent the mass of each gas and 
DM particle in the highest resolution region, and Columns 6 
and 7 are the softening lengths in that region. The cosmolog- 
ical parameters for the simulation are presented in Columns 
8-10. flo is the total matter density fraction, ho is the Hubble 
constant (100 kms^^Mpc^^) and f^b is the baryon density 
fraction in the universe. The age of the Universe is 13.5 Gyr 
in our simulation. 

To identify the main progenitor gala:xy, a friends-of- 
friends (FOF) group finder is used at regular time intervals 
in the simulation. Specifying a linking length b and identi- 
fying all pairs of particles with a separation equal to or less 



Figure 2. Star formation rate (SFR) history of disc stars . Two 
periods of star formation are visible. 



than b times the mean particle separation as friends, stel- 
lar groups are defined as sets of particles connected by one 
or more friendship relations. The other main parameter in 
the FOF algorithm is the minimum number of particles. By 
setting this parameter sufficiently high, one avoids including 
spurious objects that may arise by chance. In our simulations 
we use a linking parameter 6 — 0.01 and a threshold number 
of particles of 100. We define the largest group which has the 
highest number of the FOF identified particles as the main 
progenitor galaxy. 

Note that the ver sion of GCD+ used in this paper 
(Kawata & Gibson"2003) applies the SNe II yields calculated 
by Wooslev fc Weaver, (1995. ) . The iron yield is ambiguous 
for the SNe II nucleosynthesis model, as there is a large un- 
certainty in its value. It is well kn own that the iron yield 
shown in IWooslev fc Weaved l| 19951 ) seems to lead to lower 
[a/Fe] values, compared to those observed in low-metallicity 
stars in the solar neighbo rhood. Therefore, some au thors 
commonly use half of th e I Wooslev fc Weaver! ll 19951) iron 
yield (e.g. iTimmes etlll Il995l : iGibsonI Il997l : iGibson et al.l 
1997). The version of GCD-I- us ed in this paper app l ies th e 
actual value of the iron yield in IWooslev fc Weaver! l|l995h . 
We, however, allow for this by only comparing the relative 
difference of [a/Fe] among different samples of stars within 
our simulated galaxy. Also note that although we disperse 
metals to neighbour gas particles when a star dies, weighted 
by a kernel, metal diffusion between gas particles was not 
considered in this simulation. Therefore, the spread in the 
metallicity distributions will be artificially high, while the 
peak of the distributions should be robust. Since we likely 
overestimate the scatter in our results, we take the median 
value for approximately every 100 stellar particles in each 
[Fe/H], [O/Fe] or radius bin in Section ED and Section [33] 
respectively and plot it as a single point. This method re- 
veals more clearly any trends present in the data. 
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1 r 




Formation Time (Gyrs) 

Figure 3. [Fe/H] (upper) and [O/Fe] (lower) against formation 
time {tf) for our disc stars. Red filled and blue open circles are 
in situ and accreted disc stars respectively. 

We identify disc particles using the simulation output 
at the end of the simulation. We use the output at z = Q.l 
(t ~ 12.3 Gyr), as going to any lower redshift results in an 
unacceptable amount of contamination from low-resolution 
particles in our simulated galaxy. To define the disc, we first 
set the disc plane of our galaxy to be along the gas axis 
(x-y plane) and the rotation axis to be the z-axis. The disc 
was defined as extending radially between 4 and 10 kpc in 
the galactic plane and ±1 kpc in the z-direction. We did 
not include in our sample any stars found at less than 4 
kpc from the centre of the galaxy to minimise possible con- 
tamination with bulge stars. Furthermore, only stars with a 
rotational velocity between 150-350 kms~^ were included in 
our disc sample. This velocity range was a somewhat ar- 
bitrary choice. These combined criteria, however, worked 
well to isolate disc stars. In Fig. (TJ we show the histogram 
of rotation velocity for the stars in our disc region. Fig. [l] 
shows that the rotation component is clear and dominant in 
the disc region. Nevertheless, our sample does include some 
bulge stars, because the simulated galaxy has a large frac- 
tion of bulge stars. However, the rotation velocity criterion 
restricts the contamination only to the population of bulge 
stars which have a significant rotation velocity. Also note 
that in our si mulated galaxy the b ulge stars formed only at 
early epochs (iRahimi et al.|[201Gl ). so there should not be 
any contamination for the younger stars that are the focus 
of this paper. 



3 RESULTS 

Fig. [2] shows the total SFR as a function of time for our 
sample of disc stars. Note that Fig. [2] does not include parti- 
cles which were born in the disc region, but at this particular 
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Figure 4. Same as in Fig. [3] however now plotting only the 
median and mean every 100 particles to make the trends in our 
data clearer. The open and filled circles represent the median 
and mean for in situ stars respectively and the open and filled 
triangles the median and mean for accreted stars respectively. 




2 4 6 8 10 12 

Time (Gyrs) 

Figure 5. SFR history of our defined disc stars (upper panel) 
showing the contributions from in situ and accreted stars (red 
solid and blue dot-dashed lines respectively). In the lower panel 
we replot the SFR for our sample of disc stars but now in the 
region < -Rgc < 10 kpc and |2| < 1 kpc, with a galactic rotation 
velocity between 150 and 350 kms~^. 
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Figure 6. Face-on (upper) and edge-on (low er) view of the evolution of the galaxy from t = 7.7 - 12.1 Gyr, colour coded by the expected 
Vband luminosity llBruzual fc CharlotiliooS) . The brightest regions are the densest. The knots are due to young bright particles. 
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Figure 7. Surface density profile of stars (left) and gas (right) in the galaxy from t = 7.7 - 12.1 Gyr. 



time are not in the disc. We expect that a significant amount 
of stars, especially old stars, are kinematically heated up 
after they formed in the disc, due to our poor resolution. 
Therefore, we likely underestimate the population of disc 
stars, and our results focus on the population of relatively 
cold disc stars at the final time step. From Fig.[2]we see that 
our simulated galaxy had two episodes of star formation, the 
first occurring at an early epoch centred around 2 Gyr and 
ending before 4 Gyr. The other main episode of star forma- 
tion occurred at much later times and was more intense and 
lasted longer. The major episode of star formation respon- 
sible for making the disc started after 7 Gyr and lasted till 
the final timestep. Note that our simulated galaxy stops at 
z = 0.1. Also note that the overall star formation history is 
not like that inferred from solar neighbourhood stars in the 



Milky Way, although ther e is no consensus as to what that is 
exactly (e.g. iBertelli fc N asi 2001; Rocha-Pinto ot al. 20(51; 
lAumer fc Binnev 20091 ). This also shows that our simulated 
galaxy is different from the Milky Way. We bear this in mind 
in the following discussion. 

Our simulated galaxy formed through hierarchical clus- 
tering, and early disc formation is associated with a series 
of mergers at early epochs when bulges are also built up 
jRahimi et al.ll2010l ). One of the aims of this study is to 
show how we can infer such a merger history and building 
up history of the disc from the present-time properties of 
disc stars. To this end, first we separate the disc stars into 
accreted stars and stars formed in situ. From the current 
chemical and kinematical properties of disc stars, we study 
how chemical and kinematical properties differ between ac- 
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creted and m situ stars in Section 13.11 In Section 13.21 we 
split our in situ disc stars into two samples based on the 
age of the stars. We then compare the properties of our two 
groups of disc stars at different galactocentric radii to see 
if we can find any metallicity gradients with galactocentric 
radius and how they relate to the formation history of the 
discs. Finally, in Section |4j we summarise our findings. 

3.1 Accreted and in situ stars 

In this section, we compare the chemical and kinematical 
properties of accreted and in situ stars. We trace back the 
formation time and location for all the disc stars, and any 
stars that are born within a radius of 20 kpc from the largest 
progenitor and end up within the area defined as the disc 
at the final timestep are given the title "in situ stars" . Stars 
that form at a radius greater than 20 kpc from the cen- 
tre of the galaxy and end up within the disc at the final 
timestep, we term "accreted stars". We chose 20 kpc arbi- 
trarily because we discuss only rough trends between the 
two populations. We experimented using a larger radius up 
to 30 kpc, and generally found the same conclusions. 

Fig. [3] shows [Fe/H] and [O/Fe] versus formation 
time for accreted and in situ stars. The evolution of disc 
metallicity with t ime has been obse r ved in disc stars of 
the Milky Way (Rcddv ct all l2003l : iBensbv et al] |2004| : 
iBernkopf fc Fuhrmann.2006il although only with small num- 
ber statistics. Oxygen is one of the so-called a-elements. 
These elements are primarily produced in massive stars 
with short lifetimes that explode as SNe II. Iron is pro- 
duced predominantly in SNe la, from lower mass binary 
st ars with longer lifetim es (we apply the model proposed 
bv lKobavashi et al.ll200oh . Accreted stars have lower [Fe/H] 
and higher [a/Fe] since in this particular galaxy they form 
in early epochs before the enrichment from SNe la becomes 
important. In situ stars continue to be born up until the final 
time-step of the simulation. Fig. |3] shows the same informa- 
tion but taking the median and mean every 100 particles to 
show the trends more clearly. There is only a very small dif- 
ference between the two methods (especially for the younger 
stars which are the focus of this paper). 

In the upper panel of Fig. O we replot the SFR history 
for our sample of disc stars now including the individual 
contributions from m situ and accreted stars. We see that 
the accreted stars contribute prominently to the fraction of 
old stars in our sample. At later times we see that the disc 
grows due to the formation of stars in situ in the disc. In 
situ stars therefore are mainly responsible for the formation 
of the disc at times later than 7 Gyr. This is closely re- 
lated to the formation time of the bulge of this galaxy, which 
takes place over the first few Gyr (see Fig. 1 of lRahimi et al.l 
I2OIOI ). We can therefore deduce that the disc of our simu- 
lated galaxy only really begins to build up after the bulge is 
fully formed. However, Fig. [S] also demonstrates that some 
disc stars can form in situ during the formation epoch of 
the bulge stars due to mergers. It is possible that some of 
these stars may have formed in the vicini ty of the bulge, but 
ended up in the disc at the fina l timestep l|Roskar et al.|[200^ : 
ISanchez-Blazguez et al.ll2009l ). In the lower panel of Fig. O 
we replot the SFR history as in the upper panel, but now 
including the central 4 kpc of the bulge component. We see 
that star formation in the inner region starts earlier than in 
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Figure 8. [O/Fe] versus [Fe/H] at the final timestep. Red filled 
and blue open circles are in situ and accreted stars respectively. 
As expected from Fig. [3] accreted stars have higher [O/Fe] ratios. 



the outer region. In Fig. [6] we show how our simulated galaxy 
evolves from t = 7.7 - 12.1 Gyr. We can clearly see that be- 
fore 8 Gyr the radial extent of the galaxy is less than 4 kpc. 
Only after 8 Gyr does the disc grow beyond 4 kpc. This ex- 
plains why in Fig.[2]and the upper panel of Fig. (5] it appears 
as if the galaxy suddenly starts forming stars at around t 
= 8 Gyr. It is simply due to our selection criteria for the 
extent of the disc region. Fig. [7] shows the surface density 
radial profiles of both the stars and gas for our galaxy cor- 
responding to the times given in Fig. [6l The stellar surface 
density of the inner region {Rgc < 4 kpc) increases earlier 
(t < 8.3 Gyr) than the outer region. From these figures, we 
deduce that the simulated galaxy formed inside-out. In our 
defined disc region, the gas density is greatest at 10.3 Gyrs. 
This corresponds to the peak of the star formation rate in 
the disc. 

In Fig. |S]we plot the median of [O/Fe] versus [Fe/H], 
using the median value for every 100 particles as described 
in Section[2] Here, the abundance of oxygen is used to repre- 
sent a-elements. Accreted stars have higher [a/Fe] and lower 
[Fe/H] compared to m situ stars, corresponding to them be- 
ing formed earlier. We find that 90% of the accreted stars fell 
into the region within a radius of 10 kpc before t = 8 Gyr. 
This is before the disc starts forming at radii greater than 4 
kpc (Figs.[3]and[5}. As a result there exists a distinct popula- 
tion of accreted stars with high [g/Fe]. These ac creted stars 
likely end up in a thick disc (|Abadi et al. 2003b : .Helmi et al.l 
I2OO6I ). In fact, the velocity dispersion of the vertical compo- 
nent of the disc is 96 and 58 kms~^ for accreted and in situ 
stars respectively in the simulated galaxy. Note that this is 
significantly higher than what is observed around the solar 
neighbourhood (e.g. iHolmberg et al.ll2009t ). 

Our chemodynamical simulation demonstrates that if 
such accreted stars are formed exclusively at early epochs, 
then an accretion origin for the thick disc can explain 
the observed distinct difference in [g/Fe] between the 
thick a nd th in disc in th e Milky Way (e.g. iGratton et al.) 
1996I. I2OOOI: iFuhrmaiir] 1 19981: IProchaska et al.l I2OO0I: 



Tautvaisiene et al.' '200l'; Feltzing ct al' '2003"; 'Rcddv et alj 
2003. : . Schroder fc Pa gcl 2003; Bcnsby ct al. 2005; Fuhrman^ 
2003). We note however that it is more difficult to explain 



the thick disc stars with higher metallicity as the accreted 
stars, because the accreted stars formed in smaller galaxies 
that are likely to be more metal poor. Although our simu- 
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Table 2. Abundance gradients for various element 
abundance ratios for the intermediate and young stars 



Abundance 


Trends (dex kpc 




Abundance ratio 


Intermediate 


Young 


[Fe/H] 


-0.050 


-0.066 


[O/Fe] 


-0.005 


0.009 


[Mg/Fe] 


-0.004 


0.012 


[Si/Fc] 


-0.004 


0.006 


[N/O] 


-0.005 


-0.012 



lation shows high metalhcity and low [O/Fe] accreted stars, 
since we assume a negligible SNe feedback model (Sec. 2) 
and ignore the UV background radiation, our simulation 
likely overestimates the metallicity of the accreted stars 
formed in such sm all galaxies (e.g. iGovernato et al.1 l2007l : 
iBrooks et all l2007l ). The majority of thick disc stars in 
the Milky Way are relatively metal rich, and may require 
more intense star formation which may be more readily 
associated with the in si t u population (e.g. IChiappini et al.l 
I2OOII : iBrook et al.ll2004bl : ISchonrich fc BinnevlbOOObf T 



3.2 



Radial trend of chemical properties and age 
dependence 



In our simulation, we follow the formation history of our 
galaxy from very early times to the final timestep of the 
simulation {z = 0.1). Since the differently aged populations 
in the galaxy likely hold different memories of the formation 
history, we analyse the properties of stars with different ages. 
Fig. [2] clearly shows that our simulated galaxy has two peri- 
ods of disc star formation. The second and major episode is 
of particular interest as this is what predominately leads to 
the formation of disc stars from smooth gas accretion. Prior 
to this epoch, multip le mergers prevent smooth gas accre- 
tion to form a disc (iRahimi et al.|[2010l ). The second and 
major period lasts from t = 7 to 12 Gyr and is characterised 
by a rapid rise and fall in the SFR. It would be interesting 
to consider the differences between the stars formed at the 
two different parts of this period of star formation. Also, 
since we later would like to compare our results with the 
relatively young stars found in the Galactic disc, we there- 
fore decided to divide the stars formed in situ during the 
second episode into "intermediate" disc stars (7 < i/ < 10 
Gyr) and "young" disc stars formed during the final 2 Gyr 
(tf > 10 Gyr). Note that in this section, we only focus on 
in situ stars, to track the building up history of the disc. 

In the top panel of Fig. [9] we plot [Fe/H] against ra- 
dius from the galactic centre, Rgc, for intermediate-age and 
young stars. Here again we have plotted the median for ev- 
ery 100 particles in each radial bin. The straight and dashed 
red lines represent the best fit to the data. Here we only 
consider the disc stars within 10 kpc, because in the outer 
region there are too few particles to represent the lower den- 
sity regions of the disc. We clearly see a negative [Fe/H] 
slope with Rgc- The younger the stars, the higher the me- 
dian metallicity at any given Rgc- Table 2 shows the fitted 
slope to the median abundances for the intermediate and 
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Figure 9. [Fe/H], [O/Fe], [Mg/Fe], [Si/Fe] and [N/O] versus 
galactocentric radius. The black open and blue filled circles rep- 
resent intermediate-age and young stars respectively. The red 
dashed and solid lines are the fitted linear function to the median 
of the intermediate and young stars respectively whose slopes are 
listed in Table 2. 



young stars. Note that the exact value of these gradients 
are not important and depend slightly on how we measure 
or sample the data. Below we only discuss qualitative trends 
of these gradients. The intermediate and young stars have 
a similar slope in [Fe/H], although for the young stars the 
slope is slightly steeper. Therefore we find that during the 
major disc formation phase there is little evolution of the 
metallicity gradient, although the metallicity increases with 
time. This s eems to be consi s tent w ith what was found by 
model A of IChiappini erall l|200lh in the latter epoch of 
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Figure 10. SFR versus time at different radial regions along the disc plane 
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Figure 11. [O/Fe] versus galactocentric radius at 11.2 Gyr. The 
blue open circles and red dots represent gas and young stars re- 
spectively. The small red dots represent the yield from the disc 
stars in the last 1 Gyr. The red dashed line redraws this yield 
excluding the contribution from SNe la. 



evolution. However note that our cosmological simulations 
involve more complex processes, such as radial mixing of gas 
and stars, and a complex gas accretion history. For the young 
stars of the disc we clearly see a negative slope in [Fe/H] 
with Rgc- This kind of detailed abundance observations 
along the disc is only really available for o ur Milky Way. 
Interestingly, various observational studies l|Twarog et al.l 



Interestingly, various observational studies l | lwarog et al. 
19971: iFriel et al.ll20ql lYong et al.ll2005l: ISestito et al.ll2008l : 



Magrini et al.ll2009al : fPancino et al.ll2010l f have found simi- 
lar abundance trends for young stars to what we see in our 
simulated galaxy. 

The second row of Fig. [9] shows [O/Fe] against Rgc us- 
ing the same plotting conventions and Table 2 shows the 
slope in the median abundances for both the intermedi- 
ate and young stars. Generally speaking we see a slope 
which is close to being consistent with flat. However for the 
young stars, we do see a slight positive trend in the me- 
dian. This positive slope is also observed for Galactic open 
clusters in lYong et al] and lPancino et all (|2010l ). One 

should note that th e age range f or the youn g stars roughly 
corres ponds to the lYong et al] (|2005l ) and Pancino et al 



open clusters. Also note however that the lYong et al 



I) sample had a stronger positive gradient than the 



iPancino et al] (|201(]| ) sample, and that the actual slope is 



still uncertain. IChiappini et al] (|200ll ) showed that infalling 
gas from the halo can impact the metallicities in the outer 
regions of the disc. In Fig. [TT] we analyse [O/Fe] for gas 
and stars younger than 1 Gyr at t = 11.2 Gyr when young 
stars in Fig. [9] are forming. Fig. [TT] shows that stars inherit 
their abundance patterns from the gas. We also plot [O/Fe] 
for the yield in the last 1 Gyr from our disc stars, i.e. stars 
at l^l < 1 kpc and 4 < Rgc < 10 kpc and rotating with 
150 < Vrot < 350 kms~^ More than 70% of stars in the 
disc region, i.e. 4 < Rgc < 10 kpc and |z| < 1 kpc satisfy 
our criteria for the rotation velocity of the disc stars, yet 
93% of the mass fraction of metals is ejected from such disc 
stars, due to the young age of this population. It is therefore 
apparent that the yields from our disc stars are the domi- 
nant source of chemical enrichment, and the yield from the 
spherical component is negligible in our simulatiorQ. 

Abundance ratios, such as [O/Fe] are very sensitive to 
the star formation histories. In order to investigate how the 
star formation history changes with Rgc, and thereby bet- 
ter understand these trends, in Fig. [10] we plotted the SFR 
against time for three different -Rgc bins along the disc. 
Our three regions, going radially outwards were defined as 
extending from 4 — 6 kpc, 6 — 8 kpc and 8 — 10 kpc. The 
magnitude of the SFR decreases with increasing Rgc- By 
8-10 kpc (3rd panel from the left) there is a significant re- 
duction in the magnitude of the SFR. More interestingly, 
however, the peak of the SFR occurs at slightly later times 
as Rgc increases. Together with Figs. [5] [6] and [7] this pro- 
vides clear evidence for the inside-out formation of the disc. 
As a result, there is a higher fraction of young stars in the 
outer regions of the disc. This can explain the [O/Fe] en- 
hancement because a more significant enrichment from SNe 
II is occurring in the outer region. On the other hand, in 
the inner region, there are relatively more old stars which 
are the precursors of SNe la (producing more Fe). Fig. [TT] 
also shows the yield coming exclusively from SNe II, which 
show a constantly high [O/Fe] at all radii. This also indicates 



^ In Fig. 11, [O/Fe] of stars, gas and yields are slightly different 
from each other. However, stellar abundances come from young 
stars which formed in the last 1 Gyr, yields are only for the last 
1 Gyr, and the gas is the abundance at 11.2 Gyr. In addition, the 
plotted stellar and gas metallicities are median values at each ra- 
dial region. Stars formed only from high-density gas. Considering 
these facts, it is not surprising that they are different. 
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that the difference in the significance of SNe la enrichment 
is driving the positive [0/Fe] gradient. Therefore we con- 
clude that inside-out formation of the disc naturally leads 
to a higher fraction of young stars in the outer regions and 
a positive [a/Fe] slope in the young population of the disc. 

In Fig. [9] we also show how [Mg/Fe] and [Si/Fe] vary 
with -Rgc and Table 2 shows the value for the median fit- 
ted slopes for the intermediate and young stars. Magnesium 
and silicon are also Q-elements, and are mainly produced 
in massive stars. Therefore we observe similar trends to the 
[O/Fe] case. However, for the young stars, [Si/Fe] shows a 
slightly fiatter slope compared to [O/Fe] and [Mg/Fe] cases. 



Iwamoto et al. 


199g|) 


Gibson et al.l ( 


1997h. 



by SNe la (we adopt the yields in 
pared to SNe II as also shown in 

We also analyse [N/0] as shown in the final row of 
Fig. [9] and the fitted slopes are presented again in Table 2. 
[N/O] is expected to show an opposite trend to the [a/Fe] 
case. Accordingly, for the young stars we see a negative 
slope that has an opposite sign to the [a/Fe] case. From 
Fig. IIOI this is because at the outer radii, a greater fraction 
of young stars exist which end up as SNe II and thereby 
produce more oxygen making [N/O] lower. On the other 
hand, more intermediate mass star s (we adopt yields in 
Ivan den Hoek fc GroenewegenI 1 19971 ) died in the inner re- 
gions, producing more nitrogen, which makes [N/O] higher. 

As mentioned above, in the Milky Way, [a/Fe] in open 
clust ers also shows a tentative trend of inc reasing with ra- 
dius (|Yong et al.ll2005l : iPancino et al.ll2010l '). If the chemical 
composition of open clusters represents the properties of the 
field disc stars at similar galactocentric radii, our simulation 
demonstrates that the observed abundance trends can be 
explained by a higher fraction of young stars in the outer 
region due to an inside-out formation of the Galactic disc. 
Note that although we rely on the yields from SNe II, SNe 
la and intermediate mass stars as explained in Section 2, 
some of these values are still controversial and there are a 
variety of yields suggested by several groups. Some differ- 
ent yield sets may be able to lead to the positive [O/Fe] 
gradient without having inside-out star formation. In addi- 
tion, if somehow the infalling gas strongly contributes to the 
chemical abundance in the gas disc, and has lower [O/Fe] in 
the inner region, this can also explain the positive [O/Fe] 
gradient. Needless to say, there are more scenarios that can 
explain the positive [O/Fe] gradient. Our scenario is not a 
unique solution, but our simulation demonstrates one pos- 
sible explanation (derived with a fully cosmological simu- 
lation) which is also naturally predicted from disc galaxy 
formation in a ACDM universe. 



4 SUMMARY 

In this study we have analysed the chemistry and the dy- 
namics of the disc stars in a cosmologically simulated disc 
galaxy. The galaxy was similar in mass and size to the Milky 
Way, and contained distinct gas and stellar disc components 
iBailin et al. 2005; Connors e t a l. 2006i). 

The simulated galaxy showed two episodes of star 
formation which led to the buildup of the stellar disc. 
The first occurred at very early epochs and was mainly 
due to accreted stars being brought into the galaxy as 



a result o f mergers. These e arly mergers also built up 
the bulge l|Rahimi et al.l [201 (j ). Therefore, our simulation 
demonstrates that some stars accreted earlier can become 
disc stars, if they happen to have the right angular momen- 
tum. Our simulation also has old stars formed in situ. It is 
interesting to note that even during the mergers that are 
mainly building up the bulge some disc stars can form in 
situ and stay in the disc for a long time. Intriguingly, in 
this galaxy, a fraction of the b ulge stars have a s ignificant 
rotational velocity component (|Rahimi et al.|[201(]| ). 

The second episode of star formation starts after the 
mergers cease and continues till the final timestep. In this 
period, smooth gas accretion builds up the disc and stars 
formed in situ. We have analysed the radial trend of [Fe/H], 
[O/Fe], [Mg/Fe], [Si/Fe] and [N/O] for these m situ stars, 
especially focussing on the relatively young stars (t/ < 2 
Gyr). Our simulated galaxy shows a negative [Fe/H] gra- 
dient with Rgc- Interestingly, the Milky Way has a simi- 
lar trend as recently obser ved for open cl usters llFriel et al.l 
20021: iPancino et aLlbOld) and cepheids (iLuck et allbOolT 



Andrievskv et al.H2004l : iLemasle etal]|2008l : iPedicelli et alj 
20091 ). More interestingly, we also found a positive [a/Fe] 



gradient for the young stars in the disc. This trend is simi- 
lar to the recently observed [a/Fe] gradient for open clust ers 
in the Milky Way (|Yong et aLlbOOSl : IPancino et alll2010l ) as 
well as to predictions from some chemical evolution models 
(e.g. Magriin ct al. 2009al ). 

We find that inside-out formation of the disc can natu- 
rally produce such negative [Fe/H] and positive [a/Fe] gra- 
dients. In our simulated galaxy, the magnitude of the SFR 
declines in the outer disc. Thus we induce a more progressed 
enrichment in the inner region and negative [Fe/H] slope. We 
also find that the peak of the SFR occurs at a later epoch 
in the outer regions of our simulated galaxy, which is clear 
evidence for inside-out disc formation. As a result, the outer 
region in the disc harbours a greater fraction of young stars 
that produce more SNe II and cause higher [a/Fe]. This 
also leads to a negative slope of [N/O], because in the in- 
ner region a greater fraction of low mass stars produce more 
nitrogen. These results demonstrate that such radial gradi- 
ents of chemical abundances are sensitive to the formation 
history of the disc. 

Note that our simulated galaxy is not a late-type galaxy 
like the Milky Way, but is more like an early type disc galaxy. 
Therefore, the age distribution of the disc stars in our simu- 
lated galaxy is very different from the Milky Way disc stars. 
However we should be able to apply our simple conclusion to 
any disc that formed in situ as a result of smooth gas accre- 
tion. Therefore, we suggest that the observed positive [a/Fe] 
gradients in the Milky Way disc stars can be explained if the 
Milky Way also experienced a clear inside-out formation and 
harbours a greater fraction of young stars in the outer re- 
gion. Although this is not a unique scenario and we do not 
reject other scenarios, this is naturally expected in a ACDM 
universe. 

Recently, IWilhams etai] l|2009l ) found clear evidence for 
inside-out disc formation in M33 and claimed that there is 
a greater fraction of old stars in the inner regions of the 
disc of M33. Various authors have found that the metal- 
licity gradient in the M33 disc decr ea ses going away from 
the centre (e.g.lMagrini et al.ll2007allbl: IViironen et al.ll2007l: 
iRosolowskv fc Simonl2008l : lRubin et al.ll2008l : lMagrini et all 
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l2009bl : ICionill2009l ). More measurements of abundance ra- 
tios across t he M33 disc would be extremely interesting 
(such as: e.g. lEsteban et al]l2009[ l. Recently, there have also 
been attempts to measure the chemical properties along 
the d isc radius for disc galaxies other than our Milky Way 
(e.g. iRvder et af] l2005l: iMacArthur etahl 120091 '). Interest- 



ingly, Yoachim fc Dalcanton ( 20081 ) find a tentative detec- 



tion of more a-enhanced populations and younger luminos- 
ity weighted mean ages in the outer disc region of the disc 
galaxy FGC 1440. 

Our results are useful for comparing observations to our 
simulation where we can trace the formation history. How- 
ever, we also admit that the current chemodynamical simu- 
lation model should be improved. For example, since we do 
not allow metal mixing between particles, we likely overes- 
timate the scatter of the metallicity distribution and [a/Fe] 
at different radii. We are now working on improving our 
chemodynamical model, and new simulations will provide 
valuable information to disentangle the formation history of 
the disc galaxies from current and future observations. 
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